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1. Introduction 

Since its first appearance in the scientific literature []]], 0|, more than fifty years ago, the rapid 
(r-) neutron-capture process remains one of the most exciting and challenging questions in nuclear 
astrophysics. In particular, the theoretical quest to explain the production of r-process isotopes, 
and the astrophysical scenario where this process occurs have not yet been satisfactorily solved 
(for a general review see, for instance, Ref. []3| Q]). R-process isotopic and isobaric abundance 
distributions are typically deduced by subtracting the calculated s- and p- process contributions to 
the observed Solar System abundances. Furthermore, isotopic abundances originated in the early 
Galaxy can be directly observed in metal-poor Eu-enriched halo stars (MPEES) (i.e. [Fe/H]<— 2, 
[Ba/Eu]<— 0.7, [Eu/Fe]>+1) These combined observations reveal disparate behavior for light 
and heavy nuclei: MPEES abundance-patterns are nearly consistent from star to star and with the 
relative solar system r-process abundances for the heavier neutron-capture elements A>130 (Ba 
and above), suggesting a rather robust main r-process operating over the history of the Galaxy. Such 
a consistent picture is not seen for light neutron-capture elements in the range 39<Z<50, as the 
Solar System Eu-normalized elemental abundances in MPPES show a scattered pattern [Q, ||,[7],[8|]. 
Studies based on anti-correlation trends between elemental abundances and Eu richness at differ- 
ent metalicities provide evidence for a second different origin of isotopes below A~130 [g, [E]]. 
Measured abundances of 107 Pd and 129 I, for A<130, and 182 Hf, for A> 130— trapped in meteorites 
in the Early Solar System formation [11, [l2|] — further reinforce the idea of different origins for 



isotopes lighter and heavier than A=130. 

Different models are currently used to infer plausible sites for the nucleosynthesis of r-process- 
like neutron-rich elements. As an example, the high entropy wind model (HEW) can reproduce the 
abundance pattern over the entire mass regime using a unique core-collapse type II supernovae sce- 



nario during its different phases [13]. Conversely, calculations based on the neutrino-driven wind 
model require three different sites: low-mass supernovae from ~8-ll M Q progenitors, normal su- 
pernovae (~ 12-25 Mq), and hypernovae (~25-50 M ) [|l4|]. Besides the underlying astrophysical 
conditions, these models differ on the nuclear physics inputs used in the calculations. As an exam- 
ple, some authors claim that the calculated abundances in the region A~110, prior to the A=130 
peak— wherein r-process models tend to underestimate the productions by an order of magnitude or 
more— can be corrected by neutrino post-processing effects JIB]]. On the other hand, calculations 
made on the basis of the classical r-process have shown that the underproduction of abundances can 
be corrected under the assumption of a reduction of the N=82 shell gap far from stability [ |T6| , |l7| ]. 
Clearly, an experimental confirmation of this quenching effect would have crucial consequences in 
the search for the r-process site. In this sense, self-consistent mean-field model calculations predict 
that such N=82 shell quenching might be associated with the emergence of a harmonic-oscillator 
doubly magic nucleus 110 Zr, arising from the weakening of the energy potential surface due to 



neutron skins [18, h_9n 



It is well known that /3 -decay properties play a crucial role in r-process model calculations [16]: 
/3 -decay half-lives (Tin) of r-process waiting -point nuclei determine the pre freeze-out isobaric 
abundances and regulate the speed of the process towards heavier elements. /3 -delayed neutron 
emission probabilities (P n ) define the decay path towards stability during freeze-out, and provide 
a source of late-time neutrons. Besides astrophysical interest, these j8 -decay properties are re- 



2 



Beta decay studies of r-process nuclei at the NSCL 



J. Pereira 



lated to the /3 -strength function Sp(E), integrated over different energy ranges, and weighted by 
its strong energy dependence (relative to the Qp value) through the Fermi phase-space function 
f(E) ~ (Qp —E) 5 . Given these dependencies, 7\/ 2 and P n can be used as first probes of the struc- 
ture of p -decay daughter nuclei in regions where more detailed structure studies are prohibitive 
owing to the low production rates and the lack of spectroscopic references for identification pur- 
poses. In particular, the j3 -decay properties of nuclei in the neighborhood of 110 Zr may be exploited 



to investigate the spherical character of this potentially doubly-magic nucleus [|18[], or even a more 
exotic tetrahedral-symmetry type predicted by some authors [EO]. Complementary to these "gross" 
probes, /3-delayed /-spectroscopy has proven to be a very powerful tool to study the structure of 
nuclei in less exotic regions. Studies of the energy-level systematics as a function of neutron num- 
ber serve to infer the nuclear structure effects that dominate the physics of r-process nuclei in more 
exotic regions. 

The present paper summarizes the experimental activities carried out at the National Super- 
conducting Cyclotron Laboratory (NSCL) aimed at studying /3-decay properties of neutron-rich 
nuclei involved in the synthesis of heavy elements. The paper is divided in two sections: In the first 
part, the different experimental setups and measurement techniques are discussed. A second sec- 
tion follows, which summarizes the main results obtained in the different experimental campaigns, 
including the nuclear-structure insights inferred from the measurements. 

2. /3-decay experiments at NSCL 

Exotic neutron-rich nuclei are produced and analyzed at NSCL, Michigan State University 
(MSU). The primary beam, accelerated in the Coupled Cyclotron Facility (CCF) at few hun- 
dred MeV per nucleon, impinges onto a production target (typically Be) with average intensities 
ranging from 1.5 pnA to 20 pnA. Fragmentation reaction products are forward-emitted and sepa- 



rated in-flight with the A 1900 fragment separator [|21|], operated in its achromatic mode, using the 
Bp—AE—Bp technique [fffi]. Two plastic scintillators located at the intermediate (dispersive) focal 
plane and in the experimental area are used to measure the time-of-flight (ToF), related to the mass 
of the transmitted nuclei. 

An Al wedge is mounted at the intermediate focal plane to keep the achromatism of the A 1900; 
energy-loss experienced by nuclei passing through this degrader system provides a further filter to 
select a narrow group of elements. The selected cocktail-beam is transmitted into the experimental 
vault located downstream the A1900. Measured energy-losses (AE) of these nuclei in a silicon PIN 
detector allow to further identify these nuclei according to their atomic number Z. The particle 
identification of the separated nuclei is verified by measuring y-lines of well-known /is-isomers 
included in the cocktail-beam, in coincidence with AE and ToF signals. 

Finally, the separated and identified nuclei reach the experimental area where they are ana- 
lyzed. 

2.1 /3 -decay half-lives — NSCL Beta Counting System 

The Beta Counting System (BCS) [^3|] consists of a stack of silicon detectors, the center piece 
being a 40 x 40-pixel doubly-sided 1 mm-thick Si strip detector (DSSD), where implanted nuclei 
are detected, together with their position-correlated subsequent /3 -dec ays. Detection efficiencies 
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of /3-decays are typically better than 30%. Downstream of the DSSD, a stack of 16-strip single 
sided Si detectors (SSSD) and a 1 cm-thick Ge crystal are used to veto light particles and nuclei 
that punch-through the implantation detector, and for calorimetry of j8 particles. The output signals 
of the DSSD and SSSD detectors are coupled to dual-gain preamplifiers designed to deal with the 
wide dynamic energy range associated with implanted heavy nuclei (GeV) and emitted j3 -decays 
(keV). The whole system is enclosed in a cylindrical, thin walled Al vacuum chamber, coupled to 
the beam line. 

Implanted nuclei are unambiguously identified according to AE, ToF and total kinetic energy 
(TKE), which is determined as the sum of energy-loss signals over all silicon detectors. The time 
and pixel where an implantation event is detected are recorded. Subsequent decay events observed 
in the same or neighboring pixels within a specified time window are assumed to be correlated. The 
distribution of time differences between implantation and decay events provides the corresponding 
j3 -decay curve for a given nucleus. Multi-parameter least-square fits are calculated for these curves 
using the Batemann equations for mother, daughter and granddaughter decays, as well as back- 
ground. Half-lives of daughter and granddaughter are fixed (when known), whereas background, 
initial number of parent nuclei and their decay half-lives are left as free parameters. Alternatively, 
a maximum-likelyhood (MLH) analysis can be used in cases where Gaussian statistics are not 
realized. The MLH method has been extensively used in previous analysis where the low produc- 
tion rates of some very exotic nuclei precluded a fit analysis of their j8 -decay curves based on the 
Batemann equation [24, 26, 27]. 



2.2 f3 -delayed neutron-emission probabilities — The Neutron Emission Ratio Observer 

Neutrons emitted after the j8 decay of a nucleus are measured at NSCL with the Neutron Emis- 
sion Ratio Observer (NERO) fl28l[. The detector consists of forty-four BF3 and sixteen 3 He pro- 
portional gas counters, embedded in a polyethylene matrix with a central hole that accommodates 
the BCS chamber along its symmetry axis. /3-delayed neutrons are moderated in the polyethylene 
block and subsequently detected by the proportional gas counters. Signals from these counters 
are registered in a multi-hit TDC during a 200 /is window following j3 decay; this time interval 
was chosen on the basis of an average moderation time of about 150 /is measured for neutrons 
emitted from a Cf source. The P n value of a nucleus is deduced from the ratio of the number 
of detected /3 -delayed neutrons— corrected for background contamination and limited efficiencies 
around 40%— to the number of j3 decays. 

2.3 j3 -delayed /-spectroscopy — The Segmented Germanium Array 

Another detector widely used in j8 -decay experiments at NSCL is the Segmented Germanium 
Array (SeGA) p9[|. Sixteen detectors from SeGA can be mounted in compact geometry around the 
BCS. The y-peak detection efficiency is 7% and the energy resolutions 3.5 keV at 1 MeV. j8 decays 
from the nuclei implanted in the DSSD can populate the excited states of the daughter nuclei. De- 
exciting 7 rays are detected with SeGA whenever a trigger signal from the /3 decay is registered. 
The combination of SeGA with the BCS provides a mean to dissentangle the level ordering of 
excited states in the daughter nucleus. The j3 branching ratios derived from 7 ray intensities can 
facilitate the assignment of quantum numbers {J 71 ) to the excited states. 
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Figure 1: (Color Online) Solar System r-process isotopic abundances (dots) compared with two classical 
r-process calculations. The first one was obtained using QRPA calculated half-lives for every nucleus (green 
line); the second one was obtained after replacing the calculated half-life of just 78 Ni by the new measured 



value (red line). Figure based on Ref. [35 1 



3. Experimental results and discussion 

3.1 Half-life of the doubly-magic r-process waiting point 78 Ni 

One of the main achievements at NSCL, regarding r-process motivated experiments, was the 
measurement of the half-life of the doubly-magic 78 Ni nucleus, along with half-lives and j8 -delayed 



neutron branchings of neighboring isotopes [|24J]. The new measured half-life, 7} / 2 = 110^° ms, 



is significantly shorter than that obtained from QRPA-type calculations [[M], [31} ]. 78 Ni represents an 
important waiting point in r-process models where the neutron-captures begin at rather light masses, 



such as those assuming non standard neutron star masses [32] or those triggering the supernovae 
explosion with the collapse of an ONeMg core [|33|]. Classical r-process calculations aimed at 
reproducing light neutron-capture elemental abundances (39<Z<50) may also be directly affected 
by this waiting point. In these scenarios, the flow of matter through the N=50 shell is constrained by 
the half-life of 78 Ni, along with those of the other two already-known waiting points 79 Cu and 8() Zn, 
and has a strong impact on the formation and shape of the abundance pattern beyond the A=80 peak. 
Indeed, classical r-process calculations have shown that the replacement of the half-life of 78 Ni, 



obtained from P. Moller's QRPA model pip, by the new (more than four times shorter) measured 
value has strong impact on the synthesis of heavy r-process nuclei, with particular emphasis for 
the cosmochronometers actinides Th and U (see Fig. [[]). Much better agreement of the calculated 
half-life with the experimental value is found in the updated QRPA-model version discussed in 



Ref. [34]. of neutrons during phases where the very r-process is not yet fully developed. 
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3.2 [3 -decay half-lives and P n values of nuclei approaching the N=82 r-process path 

The half-lives of 111115 Tc, 114 - 118 Ru, 116 - 121 Rh and 119 - 124 Pd, and P„ values of llf, - 120 Rh, 



120-I22pj an( j 124 Ag, near the r-process path, were also measured at NSCL [25]. The region 
A=112— 124 is particularly interesting as it lies in the well-known abundance trough— albeit none 
of the measured nuclei are part of the r-process path itself. The measured half-lives were compared 



with QRPA calculations Q30H using ground-state deformations, and Qp and S„ predicted by both, the 



ETFSTQ model [17], which includes quenching effects of the N=82 shell in a phenomenological 



way, and the FRDM model [36]. Although neither of these inputs provided a complete agreement 
with the data, better results were obtained when ETFSTQ was used instead of FRDM. Such result 
may indicate a weakening of the N=82 shell closure for Pd isotopes with N>74, although other 
parameters entering into the calculation of the j3 -strength function should be better known. As for 
P n values, QRPA calculations using mass extrapolations and ETFSTQ ground-state deformations 
show good agreement for most of the cases, with 120 Rh being an exception. For this nucleus, a 
discrepancy of a factor 4 is found between P n values calculated based on ETFSTQ masses, and the 
experimental limit of P n ~ 5%. This discrepancy could be due to uncertainties in the quadrupole de- 
formation and <2j3 values used in the calculations. When the new P n value is included in a classical 
r-process calculations, the 120 Sn/ 119 Sn abundance ratio increases by 40%. 



3.3 Collectivity of Pd isotopes approaching the N=82 shell: /3-delayed /-spectroscopy of 

120p d 

The structure of nuclei in the vicinity of the classical r-process path are crucial for testing the 
theoretical models aimed at calculating the relevant properties of r-process nuclei experimentally 
inaccessible. These studies are particularly important in the abundance trough region of the r- 
process mass abundances. Following this motivation, an experiment aimed at investigating the 
collectivity of Pd even-even neutron-rich isotopes via E{2\) and E{A~\) /E{2\) measurements was 
performed at NSCL [[37j]. As the low-lying yrast levels are expected to rise when approaching a 
closed shell, deviations from such a trend may reveal an enhanced collective character for these 
nuclei. 

7 rays with energies of 438 keV and 618 keV were observed in coincidence with /3 parti- 
cles following decay of 120 Rh75. These transitions were assigned as the 2\ — > 0^ and 4^ — > 2f, 
respectively, in the daughter 120 Pd74. The identified E(2f) level of this nucleus, along with its 
even-even neighbor 118 Pd72, lie in near proximity of those corresponding to 108 Pd62 and 110 Pd64, 
revealing isotopic symmetry centered in 114 Pd68 (see Fig. |2[ left). Likewise, a regular trend was 
observed when comparing those nuclei with their corresponding 4 proton-particle isotones 126 Xe72 
and 128 Xe 74 . 

As shown in Fig. ^| (right), earlier IBM-2 E(2^) systematics— calculated on the basis of boson 
number established from N=82 and Z=50 closed shells— reproduce well the whole range of Pd iso- 
topes, including the measured value in 120 Pd. This result supports the notion that the collectivity of 
Pd isotopes is systematically reduced beyond the N ~66 mid-shell [[FJ, which may have important 
consequences for the calculations of the mass and /3 -decay properties of the (potential) waiting 
point 128 Pd 82 . 
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Figure 2: (Color Online) Measured E(2+) systematics of even-even Pd isotopes (circles) compared with 
(left) 4-proton particle isotones 54Xe (black squares), and (right) IBM-2 calculations (green squares). (Fig- 
ures taken from Refs. 51). 



3.4 [3 -decay properties of r-process nuclei around N«66 mid-shell 

Two other recent experiments at NSCL include the study of N«66 mid-shell nuclei below 
A=110, and neutron-rich As and Se isotopes at the N«56 sudden onset of deformation. In these 
experiments, new half-lives were measured for the first time for 88,89 As, 90 Se, 105 Y, 10 M° 7 Zr and 
lu Mo, and new P n values were determined for 103 > 104 Y, 104 Zr and 109110 Mo. In addition, upper 
limits could be established for the P n values for 105 Y, i03,i05-i07 Zr and io8,m Mo 

From the measurements of mid-shell isotopes in the vicinity of Zr it was possible to investigate 
the systematics of quadrupole deformation parameters /3 2 with neutron number. The measured T X j2 
and P n were calculated with the QRPA model using mass extrapolations and £ 2 varied over a wide 
range of deformations. As an example, Fig |3] (left) shows the calculated T\/2 of 105 Zr, as a function 
of £ 2 , compared with the new experimental data. 

The range of prolate deformations /3 2 (derived from £ 2 ), that optimized the QRPA results to 
best reproduce the measured T\/2 of 105 Zr65, is compared in Fig. || (right) with the values predicted 
by the FRDM model. Interestingly, the deformation of Zr isotopes, expected to reach saturation 
for isotopes at the mid-shell, gets its maximum value around 104 Zr64. Beyond this isotope, the 
evolution of /3 2 with neutron number deviates from the FRDM predictions, showing a more pro- 
nounced decrease. This result may indicate a shape transition region, from highly-prolate to spheri- 
cal shapes, due to a possible double shell closure at Z=40, N=70, as predicted by J. Dobaczewski et 
al. JT^]. According to these authors, the expected doubly-magic character of 110 Zr7o arises from the 
weakening of the energy potential surface due to neutron skins near the drip line, which would be 
directly related with a quenching of the N=82 shell [18, 19|] . Alternatively, other authors attribute 



the low deformations in this region to the presence of a tetrahedral-symmetric potential minimum 
at 1 Zr-jQ [20]. Further detailed /3 -delayed /-spectroscopic studies, necessary to clarify such pos- 



sibility, will have to wait for the development of new high-intensity fragmentation-beam facilities 
based on FRIB concept. 

4. Conclusions 

R-process motivated j8 -decay campaigns have been undertaken at the NSCL in recent years. 
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Figure 3: Left: QRPA-calculated T\ji for 105 Zr as a function of £2 (solid line), compared with the measured 
value (shaded area). Right: Quadrupole deformation parameter JS2 of Zr isotopes obtained from: S. Raman 
et al. [B9fl (solid dots), FRDM calculations (empty dots), QRPA-optimized half-life (shaded rectangle). 



The high-intense beams from the NSCL Coupled Cyclotron Facility, along with the high resolution 
achieved by the A 1900 separator, were crucial for producing the very exotic species near and at 
the r-process path. Furthermore, the versatility of the BCS, NERO and SeGA detection systems 
enabled the measurement of /3 -decay properties important to r-process modeling. j8 -decay half- 
lives were deduced using two independent methods: the least-square fitting of the decay-curves 
and the Maximum Likelihood technique. Delayed neutrons and 7 rays detected in coincidence 
with p particles provided a clean method to extract P n values and to analyze the low-energy spectra 
of the daughter nuclei. 

The half-life of the doubly-magic 78 Ni nucleus was measured for the first time. Classical r- 
process calculations have shown that the inclusion of the short new half-life (T 1/2 = 110+^° ms) 
of the waiting-point 78 Ni has a strong impact on the synthesis of heavy r-process nuclei, with 
particular emphasis for the cosmocronometers actinides Th and U. 

New T 1/2 of 115 Tc, 116 - 118 Ru, 121 Rh and 121 ~ 124 Pd; P n values of 118 - 119 Rh, 124 Pd and 124 Ag; 
and upper P n limits of 116 > 117 ' 120 Ru and 12 °- 122 pd were also measured. QRPA-calculated half-lives 
of Pd isotopes show better agreement with the new measured data when ground-state deformations, 
Qp and S n are taken from ETFSI-Q. Furthermore, classical r-process calculations have shown that 
the abundance ratio 120 Sn/ 119 Sn is strongly affected by the measured P„ of 120 Rh. /3-y spectroscopy 
of 120 Rh provided first information on low-energy yrast states in 120 Pd. The systematic variation 
of £"(2^) and £(4^) in the even-even Pd isotopes up to 120 Pd74 show good agreement with IBM-2 
calculations which assume good shell closure at Z=50 and N=82. 

New half-lives were also measured for 88 89 Ar, 90 Se, 105 Y, 106 107 Zr and ni Mo, along with P n 
values of 103 > 104 Y, 104 Zr and 109 - 1I0 Mo, and upper P n limits of 105 Y, 103,105-107^ and i08,m Mo 
QRPA calculations using j3 2 parameters optimized to reproduce the new measured half-life of 
105 Zr65 show a reduction of the quadrupole deformation for these nuclei. Such trend may indi- 
cate the presence of a spherically symmetric 110 Zr7o predicted by some authors. 

Extension of j3 -decay measurements to full r-process pathway requires new high-intensity 
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fragmentation-beam facilities like FRIB. In the mean time, fully microscopic j3 -decay models are 
necessary to make reliable extrapolations of /3 -decay properties to r-process nuclei. 
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